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Abstract Based on the strong coupling expansion, we 
reinvestigate the scaling behavior of the susceptibility 
X of two-dimensional O(N) sigma model on the square 
lattice by the use of Pade-Borel approximants. To ex- 
ploit the Borel transform, we express the bare coupling 
g in series expansion in x- At large N, Pade-Borel ap- 
proximants exhibit the scaling behavior at the four-loop 
level. Then, the estimation of the non-perturbative con- 
stant associated with the susceptibility is performed for 
N > 3 and the results are compared with the available 
theoretical results and Monte Carlo data. 
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1 Introduction 



Pade-Borel 
susceptibility 



In the lattice formulation of physical models, the lattice 
spacing a is a basic parameter under feasible theoretical 
control. It is therefore natural to describe observables 
in terms of a. We concern here with the possibility of 
the use of large a expansion for the approximation of 
the continuum scaling. 

In a recent paper [T], continuum scaling of lattice 
field theories was studied from the view point of strong 
coupling expansion. The model studied there is the non- 
linear O(N) sigma model on two-dimensional square 
lattice. This model enjoys properties similar to Yang- 
Mills model such as asymptotic freedom and dynam- 
ical mass generation [5] and provides us a convenient 
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testing ground of a new computational scheme. Since 
the large lattice spacing means the large bare coupling 
due to asymptotic freedom, large a expansion is equiv- 
alent with the strong couping expansion. It was then 
demonstrated that the Pade-Borel approximants of the 
large a (a: the lattice spacing) series of bare coupling 
g shows four-loop scaling for large N. The estimation 
of the non-perturbative constant , which enters into 
the relation of the correlation length £ with g such as 
i = Q exp(f^)( ^§)~^ (1 + • • •) where (3 = l/g\ 
was then estimated at N > 3. The result of estimation 
is in good agreement with the result via thermodynamic 
Bethe ansatz due to Hasenfratz et. al [3]. 

Let us briefly survey the outline of [I] . In [TJ , the ba- 
sic variable by which the bare coupling is expressed was 
the momentum mass M defined by the zero momentum 
limit of two-point correlation function. M is the practi- 
cal replacement of a. In the Fourier transformed form, 
M appears in the correlation function as 



/3(<X • <X X ) = N 



d 2 p Zexp(— ip • n) 
_ w (2tt) 2 M+p 2 + 0(p 4 )' 



(1) 



Z denotes the wave function renormalization and the 
coefficient of p 2 is rcscalcd to be one. Conventionally, 
M is considered as function of j3. In [T], however, an- 
other approach was employed to access the scaling be- 
haviors: Mutual roles of j3 and M was inverted and j3 
was expressed in M in the form j3 — Y^,k=i &kM~ k . 
Then, recalling that one must consider the behavior of 
(3 at 1/M 1, Borel transform with respect to M 
was attempted. Borel transforming the relation giving 
/3 = J2k=i( a k/k*)M~ k (bar variables denote the Borel 
transformed ones) which is entire series and improving 
the series by Pade approximants, the non-perturbative 
constant C{ has been estimated. We would note that 
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Borel transform in our use has connection with the scal- 
ing transformation, as illustrated in [Tj. 

When the second order phase transition is under 
consideration, the choice of the momentum mass as a 
basic variable describing the system seems to be natu- 
ral, since it spells the typical cooperative length of the 
system. There is, however, another variable which may 
play a similar role with 1/M or £, the susceptibility 
In this work, the susceptibility \ is defined as the sum 
of all two-point correlation functions times /3, 



(2) 



and enters into the two-point correlation function as 
d 2 p exp(— ip ■ n) 



/3(er ■ er x ) 



N 



V + o(p 4 ) ' 



(3) 



The susceptibility x describes the magnetic response 
of the system when infinitesimal external fields are ap- 
plied. Actually, from comparison of JT]) and ©, it is 
clear that the limit x ~^ 00 means the divergence of 
the dominant length scale and corresponds to the con- 
tinuum limit, thus x playing a similar role with 1/M. 
In addition, the computation of the non-perturbative 
constant C x associated with \ itself is also of physical 
interest (the definition of C x will be given in the next 
section). 

The purpose of the present paper is to apply the 
same approach taken in [T] to the susceptibility of 0(N) 
sigma model on square lattice and examine the predic- 
tive power of Pade-Borel transform on the strong cou- 
pling series for the approximation of continuum limit 
and non-perturbative quantity C x . In the present pa- 
per, we confine ourselves with N > 3. 

In the next section, we survey the series expansion 
at weak and strong couplings. Then we make attempt 
to exploit strong coupling series to access the scaling 
behavior and compute the value of C x . For the pur- 
pose we use Pade-Borel approximation method. Some 
remarks on the roles played by (3 = 1/g 2 and x are also 
given. Conclusion is given in the last section. 



2 Series expansions at weak and strong 
couplings 

The standard action of the two-dimensional lattice non- 
linear sigma model with O(N) symmetry reads 

S = -^^^„^„ + e„ (4) 

n /i=l,2 

where ei = (1,0), e2 = (0, 1) and 
1 



r 



(5) 



The vector er = (a\,(T2, ■ ■ ■ ,&n) is constrained to sat- 
isfy at every sites, er 2 = N. 

From perturbative renormalization group at N > 3, 
X near the continuum limit is found to behave as 

AttNB, {2nNf3\ -(n+i)/(n-2) 
X ~ f3C x exp(^^)^^3 -) 



f, h b 2 



(6) 



"N-2 
I + P 

where b\ and 62 represent the three- and four-loop con- 
tributions, respectively, and obtained in [HG3IMZI as 

1 



61 



N(N - 2) 
1 



-0.1888 + 0.0626A), 



z N 2 (N~2) 2 

x (0.1316 + 0.01877V - 0.02027V 2 - 0.0108A 3 ). (7) 

The multiplied constant C x is not analytically known. 
Only its large N expansion to the first order is analyt- 
ically obtained as [6] 

^-nt'-T +<*"-»>)■ w 

The leading 1/N correction to C x is big and the value 
l/(167r) in the large N limit may not become an ap- 
proximation for moderate values of N. 

At strong coupling the susceptibility is expanded in 
powers of j3. Butera and Comi [8] obtained x U P to /3 21 . 
To several orders x is written as 



72 + 32A 



X = /3(l + 4/3 + 12/3 2 + y + 2 

8(284+ 1477V + 20iV 2 ) , 
-I — : '-p? _i_ 

(iV + 2)(A + 4) P 
The result of inversion then reads 



/3 3 



200 + 6iV 
iV + 2 



P 4 



(9) 



/3 = x(l-4 X + 20 X ' 



8(29 + 147V) ^ 3 



N 



4(374+ 169A) 4 
2 + N X 
8(5212 + 3451iV+ 538iV 2 ) 



x 5 + ---). 



(2 + A)(4 + iV) x '' "° ! 

To the 21st order, the sign of coefficients of the series 
(fTO)) is alternative for all N > 0. Based upon the above 
series effective at large lattice spacing, we attempt to 
recover the scaling behavior of f3 and then estimate the 
non-perturbative constant C x for N > 3. 



3 Analysis by the use of Pade-Borel method 

In this work, Borel transform with respect to z has the 
meaning generalized to act not only on the series ^2 % 
but also on logarithms z~ const (\ogz) k 7 (k = integer), 
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(logz)'[log(logz)] m (k,l = integer), exponential func- 
tions and so on. The transformation is carried out by 
taking a certain limit in delta expansion [9,. In some 
cases, it is also useful to exploit following integral rep- 
resentation, 



B[f] 



dz &> z 



m, c 



(ii) 



I (• 2iri z z 

where C^ 1 = z denotes the Borel counter part of z. The 
contour C of the integral is along with the straight line 
parallel to the imaginary axis on the complex z plane. 
All the non-analytic ingredients such as poles and cuts 
should be in the left side of the contour if they are. 

To study the behavior of (3 near the continuum limit 
we turn to the Borel transformed counter part, B[(3] = 
/3, at small and large X where X represents the Borel 
counter part of x- The Borel transformed quantities 
have no direct physical meaning. But they have in- 
herited from physical quantities information we like to 
know, such that the constant C x , logarithmic behavior 
with the coefficient (47r) _1 and so on. 

The Borel transform of /3 at small x reads 



O M 4 20 - 
I 3 = X(l - 77\X + tttX 



8(29- 



2!' 
4(374 • 



3! 

169AT) 



(2 + 7V)4 



14iV) _ 3 

X 



(2 + iV)5! 
8(5212 + 3451iV + 538iV 2 ) 



-X 



(12) 



(2 + N)(4 + N)&. 

To compare the Pade approximants of above with the 
perturbative results, we use the following four-loop re- 
sult, 



N — 2 
4nN 
-2(N 



t + 



N -2 



log 



1 



2) 2 6 1 + 91og[k] 



1 



(N-2)H 2 
-(27 + 6(N - 



(N - 2)H 

-2(7V-2) 2 &i 
1 



2) 2 &i)log[^] 



27, r l ,< 



46 2 )(iV-2) 3 



(13) 



where 
t = log 



(14) 



The result of Borel transform of each term in (|13l) reads 

S[logx] = logx + 7£;, (15) 

C( 2 ) 

S[loglog X ] = log(logx + 1e) + 2{lQg ^ + lE) 2 

+0((logx)~ 3 ), (16) 
1 



B\ 



1 



log x log x + IE 



+ 0((logx)- d ), 



(17) 



B[ 
B 

B\ 



log log x, log(logx + 7is) 



B[ 



logx 

[ (logx) aJ 
log log X i 
(logx) 2J 
(log logx) 2 . 



logX + IE 

1 

(log X + 1e) 2 
log(logX + 7-E) , 
(logX + 7s) 2 
(log(logx + 7fi)) 2 



(logx) 2 (logX + 7B) 2 

where "/e = 0.577216 • • •. Then, we arrive at 



+ 0((logxr 3 ), (18) 
0((logx)^ 3 ), (19) 
+ 0((logx)- 3 ), (20) 
+ 0((logx)- 3 pl) 



N — 2 
AttN 



t - 



N- 2 



log 



1 



C(2)\ 
2P J 



-2(7V-2) 2 6i + 91og[|j] 
(N - 2)H 



{N-2)H 2 



- 2(N - 2) 2 b 1 + (2b\ - 46 2 )(iV - 2) s 



+ (27 + 6(7V-2) 2 6 1 )log[-i] 



27, r l n2 



where 
* = log 



2ttN > 



IE- 



(22) 



(23) 



Though the transformed series (|T2"j) in X is an entire 
series, it remains to show slow convergence to the exact 
result. To remedy the situation, we further make use of 
Pade method to extrapolate the transformed series to 
the large X region. Since (3 behaves logarithmically for 
large X) it would be natural to exploit the diagonal type 
of Pade approximants, P\ m /m\- Here, f3[ m / rn \ denotes a 
rational function of X> where both of numerator and 
denominator are polynomials of degree m. Then we will 
compare the behavior of P[ m / m ] with fi cont to the four- 
loop level. From N — 3 to 20 and at TV = oo, we have 
constructed diagonal Pade approximants at 6th, 8th, 
10th 22nd orders. The 22nd order is the highest 
order available from the work of Butera and Comi [8] . 
Singularity at positive real axis of X appears at 20th 
order for N = 3, 4, • • •, 16. In other cases, there is no 
pole at X > 0. Thus the occurrence of the singularity 
at positive real axis is presumably accidental, and we 
assume that the singularity is originally absent on the 
positive real axis at all N < 3. In any case, we show 
here the quantitative results at 22nd order where the 
pole is absent on the positive real axis. 

To see how Pade approximants recover the contin- 
uum scaling, we have plotted /Srn/n] and /3 con t for N = 
3, 4, 8 and oo ( j3 cont is drawn by using estimated value 
C app for C x .). From Fig. 1, we observe the approximate 
continuum scaling at N = 8 and N = oo already around 
logx ~ 2. On the contrary, it may be hard to say that 
the scaling is seen at TV = 3 and perhaps at N = 4. 
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Though N = 3 and 4 cases lack a clear signal of scal- 
ing, we have carried out estimation of C x at all N > 3 
by the fitting of (5 C ont to /3[ m / m i . In the fitting, recall 
that there is one unknown quantity C x which we like to 
estimate. We vary the value of C x and seek the value 
at which the curve f3 cont is tangent at a point to the 
Pade approximant. We show in Tablel the estimated 
result denoted by C app at 22nd order from N = 3 to 
N = 20. Our results are slightly smaller than those of 
Butcra and Comi [5]. The point around which the fit- 
ting is realized is \ = 24.850, 14.797, 11.903, 10.558, 
9.821, 9.644, • • •, respectively for N = 3, 4, 5, 6, 7, 8, 



Table 1 Estimated values denoted as C a pp for the non- 
perturbative constant C x . On estimation we have used Pade- 
Borel approximants of ft at order [11/11]. Cbt means the 
estimation due to Butera and Comi [5]. 



N 


Capp 


Cbt 


N 




Cbt 


3 


0.00932 




12 


0.13361 


0.134 


4 


0.03273 


0.034 


13 


0.13782 




5 


0.05611 


0.059 


14 


0.14204 


0.143 


6 


0.07526 


0.077 


15 


0.14574 




7 


0.09046 




16 


0.14899 




8 


0.10258 


0.1035 


17 


0.15187 




9 


0.11304 




18 


0.15446 




10 


0.12048 


0.1212 


19 


0.15678 




11 


0.12723 




20 


0.15886 





Comparison with Monte Carlo results is available 
for N = 3, 4 and 8. At N = 3, the two results, C x = 
0.0146(±0.0010) [TO] and C x = 0.0130(±0.0005) QI] 
were reported. At N = 4, C x = 0.0329(±0.0016) in [TO] 
and C x = 0.0383(±0.0010) in [TO]. At N = 8, C x = 
0.1037(±0.0004) in [TO] and C x = 0.1028(±0.0002) in 
[11] . Apparently our estimation at N = 3 is much 
smaller than any of Monte Carlo data. At N = 8 on 
the other hand, our result is consistent with the Monte 
Carlo results. 

One of sources of discrepancy is the cut-off effects 
at weak coupling. Formal expansion of the standard 
action yields terms irrelevant in the continuum limit 
and these terms obscure the asymptotic scaling. For 
example, according to [TO], the leading correction to 
the asymptotic scaling of perturbative renormalization 
group result (TTO1) becomes larger for smaller N. This 
may explain the large discrepancy at N = 3. Unfortu- 
nately, since concrete information is not available yet, 
we cannot continue further. In the large N limit, how- 
ever, we can discuss on the issue quantitatively. The 
discussion will be presented later. 
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Fig. 1 First three plots show four-loop perturbative j3 (dot- 
ted lines) and Pade -Borel approximants at order [11/11], 
respectively for JV = 3, 4 and 8 (solid lines). The last plot 
shows p at N = co in the continuum limit (dotted line) and 
Pade -Borel approximants at order [11/11] (solid line). 

Now, we turn to show that one can gain rough range 
of uncertainty of the estimation by the investigation of 
the near diagonal Pade-Borel approximants. As previ- 
ously stated, known logarithmic behavior of j3 at large 
enough x selects the diagonal Pade as the most appro- 
priate one. However, reliable approximants would have 
stability for small shift of degrees of numerator and de- 
nominator. Fig. 2 shows the plot of three graphs of (3 
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for N = 3 at 22nd order, /% 2 /io], /3[n/n] and /3[i /i2]- 
For /3[i2/io] and /3[io/i2]j clear sign of limitation around 
log X ~ 2 is indicated. The fitting gives C x = 0.01355 
for /9[i2/io] an( i Cx = 0.01448 for ^[10/12] ■ We have ex- 
amined the results of near diagonal Pade at 22nd order 
from N = 3 to 10. The results are summarized in Table 
2. 



0.6 
. 5 
.4 
. 3 
. 2 
. 1 





[11/11] 


[10/12] \ 


\ [12/10] 
\ N=3 



-2 2 4 6 

log X 

Fig. 2 Three graphs of j3 for N = 3 at 22nd order, £[12/10] > 
£[11/11] and £[10/12]- 



Table 2 Estimated results of C x at 22nd order for ap- 
proximants at diagonal and near diagonal assignment of 
numerator- and denominator-degrees. The last column shows 
the ratio of the average of the two off-diagonals to the diag- 
onal. 



N 


[12/10] 


[11/11] 


[10/12] 


ratio 


3 


0.01355 


0.00932 


0.01448 


1.504 


4 


0.03629 


0.03273 


0.03703 


1.120 


5 


0.05899 


0.05611 


0.05962 


1.057 


6 


0.07770 


0.07526 


0.07822 


1.036 


7 


0.09261 


0.09046 


0.09311 


1.027 


8 


0.10456 


0.10258 


0.10504 


1.022 


9 


0.11427 


0.11240 


0.11472 


1.019 


10 


0.12228 


0.12048 


0.12274 


1.017 



Note the large difference for small TV between [12/10] 
and [11/11], and [10/12] and [11/11]. For example, dif- 
ferences between [1 2/ 1 0] and [1 1 / 1 1] are . 00423, . 00356 
and 0.00198, respectively for N = 3, 4 and 8. In par- 
ticular, for N — 3, the size of C x itself is small and, 
as explicitly shown in Table 2, the ratio of the average 
of [12/10] and [10/12] to [11/11] is about 1.5. Thus, 
we find that the estimation status for N = 3 has not 
reached to enough stability yet. We have also checked 
that the proliferation range of estimated values within 
[to + 1/to — 1], [m/m] and [m — 1/m + 1] becomes nar- 
rower as the order 2to increases ( see Table 2). These 
features would lead us to conclude that, for N = 3, the 
order of strong coupling expansion is still short. 



So far, we have discussed the model at finite N. We 
now turn to the N limit. At N = 00, we find the exact 
result from (El), 



C x (oo) 



7T 

To 



0.19635- 



(24) 



Our result at 22nd order with diagonal /?[n/n] gives 
C x = 0.19970 • • • which is 1.7 percents larger than (pMf. 
One reason of the discrepancy is, of course, the trun- 
cation of large X series. However, the presence of the 
lattice artifact also prevents us from accurate estima- 
tion. As is well known, X dependence of /3 is exactly 
specified by the gap equation, 



P 



V 



1 



(2^) 2 x- 1 +2E M =i, 2 (l-cos^)- 



(25) 



On the lattice artifact, we can estimate its effect as 
follows. Near the continuum limit, we have series ex- 
pansion, 

1 



P 



4?r 



log(32 X ) 



1 



32^ X V 
(51og(32 X ) 



log(32 X ) 
-7) + .-.. 



1) 



(26) 



1024tt X 2 

The second, third and higher order terms represent the 
lattice artifacts. These are neglected in ((6]) and (fl~3l) . 
Borel transform of (l26l) reduces the effect to some extent 
as 

1 ,. ,1 5 



P 



-(log(32 x )+ 7£ ) 



(27) 



47r v-ov- /v, ■ 327] .- 10 247T X 2 

If we include the second term — or r _ we have C v = 

32-rrx * 

0.19652 by Pm/u] at X = 6.331. Third term incor- 
poration gives C x — 0.19639 at X = 5.69. Thus, much 
accurate estimation is obtained by fitting around a non- 
large value of x- 

As the last argument, let us consider what comes out 
when x is expressed in j3 as in the conventional manner. 
From the perturbative result ([B]), it follows that 



1 * 
lo g77- 



4nN (3 N + 1 
N- 2 

h 



N - 
2b 2 ~ 



lop 



2ttN 



N -2 A^-2 



log/3 



2/3 2 



(28) 



When taking j3 as the basic variable by which X is con- 
trolled, we must consider the /3 — > 00 limit. Then it 
would be nice if the strong coupling series in /3 could 
be Borel transformed and the result allows us approxi- 
mation of scaling. The right-hand side of ([2"5f seems to 
be convenient form for Borel transforming with respect 
to fi^ 1 — g 2 . However, the scheme does not work well 
compared with the inverted version which we have pre- 
sented. We like to clarify why it is better to express /? as 
a function of X than its reverse, when Borel transform 
is exploited. We focus on the large N limit since in this 
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case all of necessary things are computed exactly (Note 
that, in the N — > oo limit, X = ~tt an< ^ the following 
argument exactly applies also to M vs 23). 

First let us remind that, by taking the N — ¥ oo limit 
of ((28|) . we find that only the first two terms survives 
in ([2"8"]l. resulting log(x/C x ) -> 4tt23 - log(27r). On the 
contrary, we find from that 32x = e 47r/3 + 4(4tt23 - 
1) + 4(-13 + 12-7T/3 - 327r 2 /3 2 )e- 47r/3 + 0( e - 87l(s ) and 

fog 77" ~ 4tt/3 - log(27r) + 4(4tt/3 - l)e" 47r ^ (29) 
^x 

+4(-15 + 28tt/3 - 647r 2 /3 2 ) e - 87r/3 + 0( e - 12jr/3 ). 

The terms of exponential in (|2T)f are the cut-off effects 
which rapidly vanish in the continuum limit. The cut-off 
effects are not included in (|28p. but they are enhanced 
by Borel transformation in the following way: Suppose 
that we make the Borel transform to access the large 23 
behavior of x(/3)- Typical example of Borel transform 
may be given by B[e~ i7Tl3 ] (transform should be done 
with respect to /3 _1 ), 

^ P ] = t t ^ 1 = Jo(2^) (30) 

n=0 ^ 

where Jo denotes the Bessel function. It is also easy 
to see that 23[/3e- 47r/3 ] = ^ J Q (2y/^Kz)dz, The term 
pk e -8irp _ ^ 2) transforms in the similar manner 
and the result ends with oscillatory function. 23 [log x] at 
large 23 is also oscillatory and the amplitude is not small. 
Thus, 23 [log x] in strong coupling expansion faithfully 
recovers that oscillatory behavior and is inadequate for 
the examination of continuum scaling and estimation 
of the non-perturbative constant C x . Rather, since the 
lattice artifact is exponentially small, it is better to use 
directly the Pade method on original strong coupling 
series, for example, as performed in [8]. 



the square lattice, there exists a literature where longer 
series expansion is computed on honeycomb lattice (14| . 
We hope to report the result elsewhere in the near fu- 
ture. 
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4 Conclusion 

We have investigated the application of Pade-Borel method 
to non-linear O(N) sigma model at N > 3. By ex- 
panding the bare coupling in terms of susceptibility x 
at large lattice spacings, we have examined the scal- 
ing behavior near the continuum limit by using Pade- 
Borcl method. Scaling behavior has been observed at 
large enough N. For N = 3, we find large discrep- 
ancy between the estimated C x and its Monte Carlo 
data. For N — 4 the estimated value of C x is close 
to existing Monte Carlo data and the result in [5]. For 
larger N, scaling behavior becomes gradually clearer as 

increases and the estimated values of C x are in good 
agreement with those in literatures. 

As for small N cases, it is desirable that the longer 
series is brought to scaling examination. Though not on 



